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ABSTRACT 

Hydraulic  ram  is  the  phenomenon  associated  with  the  dam- 
age caused  by  high  speed  projectile  impacts  with  fuel  tanks. 
It  is  divided  into  five  phases,  entry,  shock,  drag,  cavity 
and  exit.   The  entry/  shock  and  drag  phases  were  of  concern 
in  this  study. 

An  analytical  solution  for  the  pressures  developed  in 
the  shock  phase  was  formulated  based  on  bullet  impact  para- 
meters. 

Pressures  developed  in  the  drag  phase  were  experimentally 
correlated  with  theoretical  predictions.   These  pressures 
were  proven  to  be  dependent,  on  front  wall  wave  reflections. 
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I.   INTRODUCTION 

The  study  of  aircraft  survivability  has  become  increas- 
ingly important  in  recent  years.   This  increase  in  impor- 
tance is  due  to  several  factors.   Among  these  are  the  esca- 
lating cost  and  concomitant  complexity  of  modern  aircraft 
making  the  loss  of  an  aircraft  less  acceptable.   The  in- 
creased complexity  of  modern  aircraft  also  tends  to  make 
them  more  vulnerable  to  anti-aircraft  weapons.   At  the  same 
time  anti-aircraft  defensive  capability  is  increasing;  wit- 
ness the  widespread  introduction  of  hand  held  infrared 
guided  missiles  of  the  Redeye  class.   These  factors  make 
survivability  analysis  and  design  of  economic  and  military 
importance.   One  of  the  most  promising  areas  for  survivabili- 
ty analysis  is  the  fuel  system;  in  particular  fuel  cells. 
This  is  because  the  fuel  cells  occupy  a  large  percentage  of 
the  total  volume  of  the  aircraft  and  are  therefore  most 
susceptible  to  damage.   Damage  to  the  fuel  cell  can  result 
in  fuel  starvation,  structural  damage,  fires  or  explosion; 
any  of  which  can  result  in  the  loss  of  an  aircraft. 

The  interaction  of  the  projectile  and  the  fluid  in  the 
fuel  cell  produces  damage  that  is  many  times  greater  than 
that  caused  by  impact  with  an  empty  cell.   One  of  the 
mechanisms  that  can  cause  damage  is  fire  and/or  explosion 
of  the  fluid.   Hydraulic  ram  refers  to  the  complicated 
dynamic  interaction  of  the  projectile,  cell  walls  and  the 
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fuel.   These  interactions  produce  internal  pressure  pulses 
which  load  the  fuel  cell  structural  members  and  produce 
damage . 

Hydraulic  ram  is  usually  divided  into  several  phases 
that  are  associated  with  different  major  pressure  pulse 
characteristics  that  occur  during  the  passage  of  the  projec- 
tile through  the  fluid.   The  five  phases  of  hydraulic  ram 
are  the  entry,  shock,  drag,  exit  and  cavity  phases.   The 
entry  and  exit  phases  occur  when  the  projectile  penetrates 
the  respective  fuel  cell  walls.   The  shock  phase  is  associ- 
ated with  the  initial  transfer  of  the  projectile  kinetic 
energy  to  the  fluid  in  the  cell.   This  energy  transfer 
creates  a  strong  shock  wave.   The  drag  phase  occurs  as  the 
projectile  continues  through  the  fluid  and  transfers  more 
energy  to  the  fluid  through  viscous  dissipation  and  pressure 
drag.   The  cavity  phase  occurs  as  the  cavity,  which  was 
created  in  the  drag  phase,  overexpands,  is  compressed,  re- 
bounds and  continues  to  oscillate  until  the  energy  is 
dissipated.   The  entry,  shock  and  drag  phases  will  always 
occur  in  order.   The  exit  phase  may  not  occur,  and  it  may 
precede  or  overlap  the  cavity  phase  depending  on  the  size  of 
the  tank,  projectile  impact  velocity  and  other  factors. 
Hydraulic  ram  research  has  been  conducted  by  various  investi- 
gators.  This  study  was  designed  to  extend  Yurkovich's 
[Ref.  1]  analytical  solution  for  the  shock  phase  fuel  cell 
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pressure  loadings  and  to  verify  experimentally  the  analyti- 
cal solution  of  drag  phase  pressures  of  Lundstrom  [Ref.  2] 
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II.   DESCRIPTION  OF  HYDRAULIC  RAM  PHENOMENON 

In  analyzing  hydraulic  ram  all  five  phases  must  be  con- 
sidered.  The  damage  associated  with  hydraulic  ram  is  depen- 
dent upon  the  anti-aircraft  threat  and  the  physical  charac- 
teristics of  the  fuel  cell.   The  threat  determines  the  path, 
energy  level  and  size  of  impacting  projectiles.   The  fuel 
cell  size,  geometry,  location  and  construction,  in  conjunc- 
tion with  the  threat,  determines  the  dominant  phase  of  the 
hydraulic  ram  phenomenon.   For  example,  there  may  not  be  an 
exit  phase  for  a  low  energy  threat  impacting  a  tank  with 
sufficiently  large  dimensions  along  the  projectile  path. 

In  any  phase  of  hydraulic  ram  the  damage  is  related  to 
the  energy  dissipated  in  that  phase.   This  makes  the  deter- 
mination of  the  energy  loss  during  the  entry  phase  (wall 
penetration)  essential.   The  energy  dissipated  during  the 
entry  phase  is  a  function  of  wall  parameters  (thickness, 
yield  strength)  and  projectile  initial  energy,  obliquity  of 
impact,  material  properties,  and  tumbling.   Studies  by  Soper 
[Ref.  3]  have  shown  that  the  percentage  of  energy  loss  for 
untumbled,  normal  entries  can  be  approximated  by: 

%  E,  ¥    =4^   =  200  gy  AT   1 
lost   E.         m      V.2 
i  l 

The  studies  also  have  shown  that  the  penetration  energy  loss 
is  independent  of  projectile  shape.   If  the  projectile  has 
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not  tumbled  and  strikes  normal  to  the  wall,  the  damage  area 
(A)  approximately  equals  the  cross  sectional  area  of  the 
projectile. 

If  the  projectile  energy  is  sufficiently  high,  it  will 
penetrate  into  the  tank  fluid  and  produce  a  strong  hemis- 
pherical shock  wave  centered  about  the  penetration  point. 
This  is  the  start  of  the  shock  phase.   The  shock  phase  is 
characterized  by  high  pressures  of  short  duration  (microsecs) 
that  are  independent  of  fuel  cell  geometry.   After  the  bullet 
enters  the  fluid,  it  has  little  effect  on  the  shock  wave 
strength  or  shape.   Investigations  have  shown  that  shock 
phase  pressures  generated  are  a  function  both  of  the  energy 
transferred  to  the  fluid  at  impact  and  of  the  fluid  proper- 
ties.  The  shock  wave  pressures  are  attenuated  rapidly  as 
the  wave  is  propagated  through  the  fluid.   Therefore,  the 
damage  associated  with  the  shock  phase  is  generally  limited 
to  the  entry  wall.   This  damage  is  caused  by  the  impulsive- 
ly generated  pressure  field  combined  with  whatever  damage 
(cracks,  etc.)  the  initial  penetration  of  the  wall  may  have 
caused. 

An  analytical  formulation  of  the  shock  phase  problem  has 
been  postulated  by  Yurkovich  [Ref .  1] .   In  this  formulation 
the  shock  formation  process  of  the  shock  phase  has  been 
modeled  as  a  point  energy  source  release.   The  formulation 
assumes  that  the  fuel  cell  walls  are  rigid,  and  that  the 
shock  radius  is  proportional  to  a  constant  power  of  time. 
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The  assumption  that  the  radius  is  proportional  to  a  constant 
power  of  time  has  been  supported  by  several  investigators 
[Ref.  4,  5,  6]. 

The  Rankine-Hugoniot  equations  are  used  to  formulate  the 
flow  properties  at  the  shock  front.   A  fourth  equation  is 
provided  by  the  Tait  equation  of  state  for  isothermal  com- 
pressible liquid.   The  fluid  property  changes  at  the  shock 
front  are  assumed  to  be  adiabatic.   If  the  shock  radius  is 
known,  a  shock  Mach  number  can  be  computed  from  energy  con- 
siderations.  From  the  shock  Mach  number  the  shock  front 
pressure  and  density  can  then  be  computed,  using  the  above 
equations.   Reference  2  postulates  the  initial  radius  is  a 
function  of  energy  deposited  given  by: 

■p 

Rn  =   (  % )1/3  (II-2) 

Po  C2  2* 

where  E   is  the  energy  released  in  the  shock  phase,  C  is  the 
speed  of  sound  and  P   is  the  density  in  the  undisturbed 
fluid. 

The  flow  parameters  behind  the  shock  front  are  com- 
puted using  the  unsteady  motion,  one-dimensional  equations, 
for  conservation  of  mass,  energy  and  momentum,  again  coupled 
with  the  state  equation.   To  use  these  equations  a  power  law 
density  distribution  was  assumed.   The  radial  velocity 
distribution  was  found  by  integrating  the  continuity  equa- 
tion.  Integration  of  the  momentum  equation  then  yields  the 
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pressure.   Integration  of  the  energy  density  behind  the  shock 
front,  knowing  the  shock  Mach  number  and  radius,  would  then 
yield  the  energy  released  during  the  shock  phase.   Since 
the  energy  released  is  assumed  to  be  known,  solving  the 
energy  integral  for  the  shock  Mach  number  is  possible.   Once 
the  shock  Mach  number  is  known  as  a  function  of  radius,  the 
shock  radius  as  a  function  of  time  can  be  found  by  integra- 
tion.  It  is  then  possible  to  determine  the  pressure,  shock 
Mach  number  and  radius,  as  a  function  of  time,  anywhere  in 
the  field.   In  this  formulation,  then,  the  energy  released 
in  the  shock  phase  must  be  known. 

Since,  in  most  cases,  not  all  of  the  projectile  energy 
has  been  released  during  the  shock  phase,  another  phenomenon 
of  hydraulic  ram  becomes  important.   This  phase  is  referred 
to  as  the  drag  phase.   It  is  distinguished  from  the  shock 
phase  by  the  gradual  acceleration  of  the  fluid  vice  the 
impulsive  acceleration  characteristic  of  the  shock  phase. 
The  energy  transfer,  therefore,  is  due  to  viscous  dissipa- 
tion and  form  drag.   Drag  is  normally  characterized  by 
defining  a  drag  coefficient.   The  drag  coefficient  of  a 
projectile  is  not  a  constant  but  varies  with  velocity, 
shape  and  tumbling  behavior.   If  the  projectile  velocity  is 
below  the  drag  divergence  Mach  number,  however,  changes  in 
velocity  have  only  a  small  effect  on  drag  coefficient  size. 
The  effect  of  shape  and  tumbling  is  then  the  dominant 
variable  for  estimating  the  drag  coefficient  for  a 
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projectile.   Experiments  [Ref.  2]  have  indicated  that  soon 
after  entry  (within  one  foot)  the  projectiles  will  begin  to 
tumble.   Furthermore,  because  of  the  high  stagnation 
pressures  on  the  nose  of  the  bullet,  soft  nosed  projectiles 
may  begin  to  deform.   This  phenomenon  effectively  changes 
the  shape  of  the  projectile  and  therefore  the  drag  coeffi- 
cient. 

The  pressure  disturbance  in  the  drag  phase  is  propagated 
through  the  fluid  at  the  fluid  sonic  velocity.   Since  the 
energy  transfer  is  more  gradual,  the  peak  pressures  are 
lower  during  this  phase.   But  the  pulses  are  of  much  longer 
duration  (hundreds  of  microsecs) .   During  this  phase  a  cavi- 
ty is  formed  behind  the  projectile.   The  cavity  is  cylindri- 
cal in  shape  and  its  size  and  expansion  rate  are  dependent 
on  the  rate  of  energy  dissipation  by  the  projectile. 

Lundstrom  [Ref.  2]    has  formulated  an  analytical  solution 
for  the  drag  phase  pressures.   In  this  model  a  velocity  de- 
cay factor  B  is  defined  as  a  function  of  the  drag  coefficient, 
mass,  area,  and  density  of  the  fluid: 

B  =  p  C   A/2m  (II-3) 

Then  the  rate  of  energy  dissipation  with  distance  is  given 
as: 

^  =  m  BV2  (II-4) 

QX 

where  B  is  assumed  to  be  independent  of  velocity.   The  flow 
field  properties  are  then  solved  for  by  using  a  potential 
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function  that  satisfies  the  wave  equation.   A  line  of  sources 
is  distributed  along  the  projectile  path  with  the  strength 
of  the  sources  being  determined  by  an  energy  balance.   The 
effect  of  wall  restraint  on  cavity  growth  is  not  considered, 
so  that  the  analysis  will  not  yield  useful  pressures  during 
cavity  collapse. 

The  pressure  wave  reflections  at  the  fuel  cell  boundaries 
are  assumed  to  be  similar  to  those  of  a  free  surface.   A  sys- 
tem of  negative  image  sources  is  utilized  which  result  in 
zero  pressure  perturbation  at  the  walls.   Pressure  waves  re- 
flected from  the  cavity  surfaces  are  not  considered. 

The  next  phase  to  be  considered  is  the  cavity  oscilla- 
tion phase.   This  is  a  complex  phase  and  the  least  under- 
stood.  Part  of  the  complexity  is  due  to  the  many  reflected 
pressure  waves  that  have  had  time  to  be  of  significance. 
The  cavity  is  formed  during  the  drag  phase  as  the  momentum 
of  the  fluid  becomes  large  enough  to  cause  the  fluid  to 
separate  from  the  projectile  surface.   A  void  is  created 
which  is  filled  with  fluid  vapors  and  air  that  enters  through 
the  wall  penetration  point.   The  cavity  then  grows  until  it 
is  arrested  by  fluid  pressures  developed  by  the  walls,  and 
collapses.   During  the  collapse  the  momentum  of  the  fluid 
causes  the  entrapped  gases  and  air  to  be  compressed  and  the 
cavity  begins  to  expand  again.   It  is  this  alternating  ex- 
pansion and  contraction  that  produces  low  amplitude  pressure 
pulses  of  long  duration.   The  phase  continues  until  all  the 
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energy  is  absorbed  by  the  fluid.   High  speed  motion  pictures 
[Ref.  7]  have  shown  that  the  pressure  oscillations  can  lead 
to  eventual  catastrophic  failure  of  walls  already  weakened 
by  the  drag  and  exit  phases. 

The  exit  phase  will  occur  if  the  energy  transferred  from 
the  projectile  during  the  entry,  shock  and  drag  phases  has 
not  dissipated  all  the  kinetic  energy  of  the  projectile  by 
the  time  the  exit  wall  is  reached.   Pressure  loadings  due  to 
projectile  proximity  to  the  wall,  and  fluid  motion  produced 
in  the  earlier  phases,  can  cause  catastrophic  damage  to  the 
exit  wall.   If  this  occurs,  the  latter  phase  (cavity  oscilla- 
tion) causes  little  further  damage  but  does  pump  fuel  from 
the  tank,  causing  fire  hazard. 

From  the  above  discussion  it  is  apparent  that  each  of 
the  different  phases  of  hydraulic  ram  produces  damage  to  the 
fuel  tank  in  a  distinct  manner.   The  entry  phase  results  in 
a  hole  in  the  tank  entry  wall.   The  damage  from  the  shock 
phase  is  usually  limited  to  the  entry  wall  and  can  lead  to 
catastrophic  failure.   The  drag  phase  results  in  lower 
pressure  amplitudes,  but  of  longer  duration,  which  can  cause 
high  levels  of  stress  on  all  fuel  cell  walls.   The  cavity 
phase  has  even  lower  pressure  amplitude  but  longer  pulse 
widths  and  is  oscillating  in  nature.   This  can  result  in 
damage  to  the  already  weakened  fuel  cell  walls.   The  exit 
phase  causes  damage  to  the  exit  wall,  usually  before  the 
cavity  phase.   It  should  also  be  apparent  that  impact  energy 
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and  the  physical  characteristics  of  the  fuel  cell  determine 
which  phase  is  dominant. 
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III.   HYDRAULIC  RAM  PRESSURE  EXPERIMENT 

The  Naval  Postgraduate  School  Ballistic  Range  was  used 
to  conduct  the  experimental  part  of  this  study.   The  basic 
components  of  the  range  were  a  rifle,  velocity  measuring 
equipment,  simulated  fuel  tank,  and  pressure  measuring  and 
recording  equipment.   Figure  III-l  is  a  schematic  of  the 
basic  components  and  Figure  III-2  is  a  down  range  view. 

The  rifle  was  a  .222  caliber  Remington.   Cartridges  were 
loaded  to  obtain  the  7,493  in. -lb.  energy  level  v/ith  a  45 
grain  spitzer  ogive  nose  bullet.   The  velocity  of  the  bullet 
was  measured  by  using  the  first  three  Avtron  No.  A914T33 
chronograph  screens  to  trigger  two  Monsanto  101B  counters. 
Two  counters  were  used  to  increase  the  accuracy  of  velocity 
measurements  and  for  redundancy. 

The  simulated  fuel  tank  was  constructed  of  welded  .160- 
inch  707  5-T6  aluminum  plate  on  the  sides,  back  and  bottom 
walls.   The  front  wall  was  .05-inch  7  07  5-T6  aluminum.   A  one- 
inch  hole  was  drilled  in  the  front  wall  for  bullet  entry. 
The  hole  allowed  calculation  of  precise  impact  energy.   The 
top  of  the  tank  was  left  open.   The  side  and  bottom  walls  of 
the  tank  were  welded  to  a  frame  of  .5-inch  thick,  three-inch 
aluminum  angle  sections  for  added  strength.   The  entry  and 
exit  wall  were  bolted  to  the  frame.   The  inside  dimension  of 
the  tank  was  18xl8xl8-inches .   Figure  III-3  shows  the  tank. 

A  Kistler  603A  quartz  pressure  transducer  was  used  to 
measure  the  hydraulic  ram  pressure  pulse  at  several  points 
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internal  to  the  tank.   The  transducer  was  mounted  in  a 
hollow  1.0-inch  stainless  steel  tube.   The  transducer  was 
held  in  the  tube  by  a  nylon  ring.   The  tube  was  then  sealed 
with  RTV.   In  this  manner  the  transducer  and  cable  were  kept 
dry.   The  tube  had  a  45   angle  at  the  tip  and  was  mounted  on 
the  tank  at  50   from  the  vertical  to  obtain  pressure  measure- 
ments along  a  radial  line.   Figure  III-4  shows  the  transducer 
and  probe  assembly.   The  probe  was  then  secured  to  the  top 
of  the  tank  so  that  it  could  be  positioned  for  measurements 
at  2,4,6,8  and  10  inches  along  a  radial  line  from  the  impact 
point. 

The  charge  signal  from  the  pressure  transducer  was  fed 
into  a  Kistler  Model  504E  charge  amplifier  with  a  545A14 
filter  and  displayed  on  an  oscilloscope.   A  Polaroid  camera 
mounted  on  the  oscilloscope  recorded  the  pressure  trace. 
The  oscilloscope  was  triggered  by  the  fourth  chronograph 
screen  mounted  one  inch  in  front  of  the  tank  wall. 

Several  firings  were  made  at  each  distance  for  consis- 
tency and  accuracy.  Typical  results  for  each  distance  are 
shown  in  Figures  III-5,  6,  7. 
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Figure  III-2.   Ballistic  Range 
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Figure  III-3.   Fuel  Cell  and  Probe  Installation, 


25 


26 


Radius  =  2  in 
Scale:   1000  psi/cmf  50  microsecs/cm 


Radius  =  4  in 
Scale:   500  psi/cm,  50  microsecs/cm 


Figure  III-5.   Pressure  vs.  Time,  En  =  7,493  in-lb  Light 


Wall. 
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Radius  =  6  in 


Radius  =  8  in 


Figure  III-6. 


Pressure  vs.  Time,  EQ  =  7,493  in- lb  Light 
Wall.   Scale:   500  PSI/cm,  50  microsecs/cm, 
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Radius  =  10  in 
Scale:   500  psi/cm,  50  microsecs/cm. 


Fiqure  III-7.   Pressure  vs.  Time,  En  =  7,493  in-lb  Light 


Wall. 
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IV.   ANALYSIS 

The  first  segment  of  this  study  was  to  improve  the  com- 
puter code  for  shock  phase  pressures  developed  by  Kappel 
[Ref .  4] .   This  code  was  based  on  the  formulation  of 
Yurkovich  discussed  in  Section  II. 

To  account  for  the  energy  dissipated  by  the  projectile 
in  penetrating  the  entry  wall  the  formulation  of  Soper  was 
added  to  the  computer  code. 

A  method  for  computing  the  energy  deposited  during  the 
shock  phase  was  needed  as  previous  formulations  assumed 
this  value  was  known.   Studies  [Ref.  4,  5]  have  indicated 
that  the  energy  released  in  the  shock  phase  is  a  constant 
percentage  of  the  impact  energy.   This  percentage  is  a  func- 
tion of  the  impact  Mach  number  with  respect  to  the  fluid 
speed  of  sound.   Plots  of  experimentally  observed  changes 
in  projectile  velocity  during  the  characteristic  shock 
formation  with  projectile  impact  were  made  by  Power  [Ref.  8] 
This  plot  is  reproduced  as  Figure  IV-i.   In  addition,  the 
corresponding  shock  phase  energy  dissipation  percentage  of 
the  impact  energy  is  shown.   It  can  be  seen  that  impact  Mach 
number  with  respect  to  the  tank  fluid  is  critical  in  deter- 
mining the  shock  phase  energy  input.   For  Mach  numbers  less 
than  sonic  the  input  percentage  is  approximately  10  percent 
and  for  moderate  supersonic  speeds  the  energy  input  percent- 
age is  greater  than  80  percent.   This  provided  an  empirical 
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procedure  for  determining  shock  phase  energy  input  as  a  func- 
tion of  bullet  impact  parameters. 

The  shock  front  pressures  are  quickly  attenuated  and  be- 
come acoustic.   Figure  IV-2  shows  the  radius  and  time  for 
the  shock  wave  to  become  acoustic  versus  input  energy  for 
both  water  and  fuel.   To  calculate  pressures  after  the 
acoustic  speed  was  reached,  a  simple  expression  based  on 
acoustic  wave  theory  was  used.   The  shape  of  the  wave  was 
assumed  to  remain  constant  and  the  pressure  amplitude  varied 
as  the  inverse  of  the  radius. 

Since  the  Yurkovich  theory  postulates  the  energy  release 
as  a  point  source,  the  model  can  be  extended  to  an  explosion 
in  the  fuel  cell.   This  type  of  event  is  characteristic  of 
an  H.E.  round  impact.   This  was  accomplished  by  integrating 
the  energy  equation  over  a  spherical  volume.   In  this  case, 
however,  the  shock  phase  energy  that  is  released  by  the  explo- 
sion must  be  known. 

Figures  IV-3,  4  and  5  show  the  predicted  pressure  dis- 
tribution for  .50  caliber,  12.7mm  and  14.5mm  AP  projectiles 
impacting  JP-5  fuel.   The  plots  illustrate  the  very  high 
pressures  generated  and  their  rapid  decay.   Figures  1V-6 
through  11  show  the  predicted  shock  phase  pressure  versus 
time  at  a  constant  radius  for  a  .222  caliber  projectile  im- 
pacting water.   Also  shown  in  these  figures  are  measured 
pressures  from  reference  9  and  analytical  drag  phase 
pressures  to  be  discussed.   From  these  figures  it  can  be 
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seen  that  the  shock  phase  pressures  could  not  be  measured. 
This  is  due  to  the  very  rapid  rise  time  of  the  shock  phase 
pressures.   The  pressure  transducers  and  associated  measuring 
equipment  were  unable  to  follow  the  signal.   This  phenomenon 
was  also  reported  by  Williams  [Ref .  10] . 

Experiments  were  undertaken  to  verify  the  drag  phase 
pressures  predicted  by  Lundstrom.   Holm  in  reference  9  re- 
ported that  the  drag  phase  pressures  could  not  be  verified. 
In  that  experiment  the  front  wall  of  the  fuel  tank  was  a 
heavy,  .5-inch,  mild  steel  plate.   Lundstrom' s  formulation 
of  the  problem  treats  the  walls  as  free  surfaces.   This  is 
done  by  using  a  set  of  negative  image  sources  so  that  the 
pressure  at  the  wall  is  zero. 

The  degree  to  which  a  wall  may  be  considered  rigid  or 
free  is  characterized  by  the  expression  [Ref.  2]: 

a  =  p  C  9/m  CIV-1) 

where  0  is  the  time  constant.   A  rigid  surface  corresponds 
to  a=0,  and  a  free  surface  to  a=».   For  a  .5-inch  thick 
steel  plate  ct=1.5,  while  for  a  .05-inch  thick  aluminum  plate 

=44.5.   The  heavy  steel  plate  is  then  an  approximation  to 
a  rigid  surface  while  the  .05-inch  aluminum  plate  is  an 
approximation  to  a  free  surface. 

A  set  of  image  sources  that  can  approximate  the  rigid 
surface  would  be  a  positive  set.   Using  this  set  of  positive 
images  on  the  front  wall  for  the  pressures  reported  in 
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reference  9  resulted  in  very  close  agreement.   This  is 
illustrated  in  Figures  IV-6  through  11  for  two  energy  levels 
The  peak  pressure  versus  radius  is  depicted  in  Figure  IV-12. 
Close  agreement  is  shown  at  medium  and  larger  radii.   The 
theory  does  underestimate  the  pressures  at  small  radii  how- 
ever. 

The  experiment  was  repeated  using  a  .05-inch  aluminum 
front  wall  to  approximate  a  free  surface.   Plots  of  the 
pressure  distribution  versus  time  at  a  constant  radius  are 
shown  in  Figures  IV-13  and  14.   These  figures  also  show 
close  agreement  between  the  experiment  and  theory  when  a 
free  surface  is  assumed. 

A  plot  of  peak  pressure  versus  radius  compared  with 
theoretical  values  is  shown  in  Figure  IV-15.   Again,  the 
theory  underestimates  the  peak  pressure  for  small  radii. 

All  of  the  theoretical  pressures  were  computed  using  a 
drag  coefficient  of  one.   This  coefficient  was  determined 
from  the  experimental  data  of  reference  3 .   The  theoretical 
formulation  allows  a  variation  of  drag  coefficient  with 
time.   This  would  account  for  any  tumbling  of  the  projectile 
No  measurements  of  the  tumbling  behavior  were  possible,  so 
a  constant  C  was  assumed.   Photographs  of  the  bullet 
trajectory  would  be  required  to  determine  tumbling  distances 

The  pulse  shape  discrepancy  was  due  to  the  assumption 
that  C  was  a  constant.   If  detailed  measurements  are  made, 
the  pulse  shape  will  be  more  accurately  predicted. 
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Reference  11  reports  that  a  data  base  was  established  for  the 
tumbling  behavior  of  14.5mm  and  12.7mm  projectiles.   From 
these  data  good  correlation  of  pressures  was  possible. 
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A  45  grain  .222  Caliber  Projectile 

□  55  grain  .222  Caliber  Projectile 

A  5.56  mm  Tungsten  Carbide  Sphere 

O  5.56  mm  Aluminum  Sphere 

V  3.17  5  mm  Aluminum  Sphere 

<]  1.59  mm  Steel  Sphere 
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Figure  IV- 1.   Shock  Phase  Energy  Release. 
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V.   CONCLUSIONS 

The  first  purpose  of  this  study  was  to  extend  the  com- 
puter code  solution  for  shock  phase  pressures.   This  was  ac- 
complished by  adding  solutions  for  the  entry  phase  and  for 
the  shock  phase  energy  release.   The  attached  code  computes 
the  shock  phase  pressures  based  on  a  knowledge  of  bullet  im- 
pact parameters  only.   It  will  also  compute  the  shock  phase 
pressures  for  a  point  explosion  in  the  fluid. 

Attempts  to  verify  the  shock  phase  pressures  experiment- 
ally were  unsuccessful.   This  was  due  to  the  inadequate  re- 
sponse time  of  the  pressure  measuring  equipment.   Until 
pressure  transducers  with  extremely  rapid  rise  times  are 
available,  shock  phase  pressures  cannot  be  experimentally 
determined. 

The  second  part  of  this  study  experimentally  confirmed 
the  analytical  solution  for  the  pressures  produced  in  the 
drag  phase,  postulated  by  Lundstrom.   In  accomplishing  this, 
the  importance  of  entry  wall  reflections  on  tank  internal 
pressures  was  confirmed.   By  using  both  a  heavy  and  light 
wall  in  the  experiment,  and  matching  these  with  appropriate 
image  sources  in  the  analytical  solution,  the  pressures  were 
correlated. 

The  use  of  a  constant  drag  coefficient  in  the  analytical 
solution  does  distort  the  pressure  wave  shape.   This  is 
because  the  use  of  a  constant  drag  coefficient  is  not  con- 
sistent with  the  tumbling  behavior  of  the  bullet.   At 
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smaller  radii  the  experimental  pressures  are  higher  than 
predicted.   This  is  due  to  the  interaction  of  the  shock 
phase  pressures  with  the  drag  phase.   That  is,  the  pressure 
transducers  are  partially  responding  to  the  impulsive 
pressure  of  the  shock  phase. 

The  results  of  this  study  show  that  the  first  three 
phases  of  hydraulic  ram  are  clearly  understood  and  can  be 
predicted  analytically.   This  should  be  of  some  assistance 
to  those  responsible  for  increasing  the  survivability  of 
aircraft  subject  to  light  or  moderate  damage. 
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C3MPUTER    PROGRAM 

c  *  * 

c  *  * 

C  *  THIS    P^OSRAM    CALCULATES    THE    FOLLOWING  * 

C  *                 PART     I  * 

C  *                          1.     SHOCK    RADIJS    VERSUS    TIME  * 

C ■  *                           2.     SHOC<    MACH    NUMBER    VERSUS    TIME  * 

C  *                 PART    II  * 

C  *                           1.     PRESSURE    VERSUS    TIME    AND    RADIUS  * 

C  *                                   FROM     IMPACT     BEHIND    THE     SHOCK  * 

C  *                                   FOR    A    GIVEN    FLJID    AND     IMPACT    ENERGY.  * 

C  *  THIS    PROGRAM    ASSUMES:     SHOCK    RADIUS    IS    PROPORTIONAL  TO* 

C  *  TIME    TO    THE    0.9    POWER    UNTIL     ACOUSTIC ,     A     POWER    LAW  * 

C  *  DENSITY     PROFILE    BEHIND    THE     SHOCK,     A    STRONG    SHOCK,  * 

C  *  ADIABATIC    ACROSS    THE     SHOCK,     THE    UNDISTURBED    FLUID  * 

C  *  PRESSURE     IS    ESSENTIALLY    ZERO,     SPECIFIC     INTERNAL  * 

C  *  ENERGY    IN    UNDISTURBED    FLUID    IS    ESSENTIALLY     ZERO.  * 

c  ^  ^ 

c  ^  ^ 

C.'-  •>/  .'-  vV  .«,   ^  ^  si,   ^*  .1,   *>.  J,   J.  ,'.   J,   J^    J.  J,   J.  ^  gu  JU    —   -'-    J-   .',   w'.  sL.    .«-  -i-   A    .',    ,'/  .'.    s1-   J^  O,  •'>   s»-  .i,   *'-   X    »•-   J,  tif  U*  •><  .V   J.  •>,  X    V-    V-  *V  «A-  -A-  J* 
*ik  'i*  y  'i'  i'  'i'  t  V   *!k  -r  t  i*  ■■■  r  'i*  t   'i-  v  v    •■•  'i*  "  -  -  ix  •.•   -  T  -t»  r  v  t  'i*  v  r  'i'  '^  'i*  -r  '<•  y  r  v  r  i*  'i'  -i'  'c  'i*  'r  v  r  t  'i'  t  t    -^  v 
C 

c 
c 
c 

C  SYMBOL    DESIGNATION 

C  WB=BULLET    WEIGHT     (GRAINS) 

C  C0=    SPEE3    3-     SOUND    FOR    FLUID     (EPS) 

C  RH00=    DENSITY     IN    UNDISTURBED    FLUID     (LB/FT3) 

C  VI=BULLET     INITIAL     VELOCITY     (FPS) 

C  DC=     INCREMENT    OF    NuN- DI MENS  I ONAL    RADIUS,     DELTA 

C  ZHETA 

C  DELT  =     INCREMENT    OF    TIME     (SECS) 

C  SIGMA=    YIELD     STRESS    OF    WALL     (PSI) 

C  R3UL=BULLET    RADIJS(IN) 

C  THIK=ENTRY    WALL    THICKNESS     (IN) 

C  XN= EXPONENT     IN    STATE    EQUATION     (7    FOR    WATER,     10.6 

C  FOR    FUEL) 

C  IFLAG=ZERO    FOR     A    HEMISPHERICAL    CASE.        IFLAG=ONE 

C  FOR    a    SPHERICAL    CASE. 

C  TMAX=MAX    TIME    FOR    CALCULATION    OF    SHOCK    PRESSURES 

C  (MIRCO    SECONDS) 

C  E03=EMERGY    RELEASE    FOR    SPHERICAL    EXPLOSIN    (INLBS) 

C  EOi=IMPACT    ENERGY     (IN-LB) 

C  EO=ENERGY    GIVEN    JP     IN    SHOCK    PHASE     (IN-LB) 

C  X=     PERCENT    EC    LOST    AT    ENTRY     WALL 

C  XX=PERCENT    3F     IMPACT    EsiERGY    LOST     IN    SHOCK    PHASE 

C  VS=     SHOCK    MACH    NO. 

C  ET=    ENERGY/VJL.     =    2-E3/(R**3) 

C  Fl=     NOiN-DIMENSIDNAL     PRESSURE    AT    THE     SHOCK    FRONT. 

C  SHOCK    MACH    NO. 

C  DODM=    DERIVATIVE    OF    NON-DIMENSIONAL     PARTICLE 

C  VELOCITY     W.R.T.     SHOCK    MACH    NO.     V S    SHOCK 

C  MACH    NO. 

C  XXM=INITIAL    MACH    NUMBER    OF     BULLET    WITH    RESPECT    TO 

C  FLUID    (V  I/CO  ) 

DIMENSION  PSI 1 (20) , Fl( 20) ,VS( 23) ,    D0DM(20) 

1 v  RXMS ( 20 ) , DR  S ( 2 0  J  , ARM  S ( 20 )  , D T ( 2  0 )  ,  T  (  2  0  )  , R  E ( 20 ) , 
1RV( 20) ,RDOD( 20) , XMSl 20 ) ,RS(20) ,RPSI( ?0) ,RF( 20) 

DIMENSION       TT(300)  ,XM(333) ,     R{ 300 ) , F2 ( 3 00 ) ,  PS  I  2 ( 300 ) , 
1D0DM2C300) ,3(300)  ,C(100)  ,F(2,20)  ,P(3,20) 

DIMENSI ON    CONRS( 303? , SHOCV( 300 ) , TOFOR( 300) 

DIMENSION    FOR, 20) ,XX( 100) ,XXM( 100) 

COMMON    BQ 

DO    500    J=l ,20 

FOR( J)=O.U 
500    CONTINUE 

N  =  19 

READ    (5,100)     W3,C0,RH00,VI , DC , DE LT ,S I 3MA ,RBUL ,THIK,XN 

READ    (5,203)     IFL AG, TMAX , E03 
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READ  (5,101) (VS( I )  T  1=1 tN  ) 

READ(5, 101  J (XX(  I  )  ,  1=1  ,32) 

READ( 5,  101 ) <XXM(  I  I  ,  1=1  ,32) 

WRITE  (6,200)  « 3, CO, RHOO, VI , DC , DELT, SI 3MA ,R BUL ,TH I K, XN 

WRITE  (6,204)  ir-LAG,TMAX 

DO  15  1=1, N 

15  WRITE(6,202)  VS(  I  ) 
DC  16  1=1,32 

16  WRITE(6,201)  XX(I),XXM(I) 
C 

C      PART  I 

C      PERCENTAGE  OF  INITIAL  ENERGY  DEPOSITED  AND  INITIAL 

C      SHOCK  RADIUS 

C 

E0I=WB*1.4286*12.*VI*VI/(64.4*10.**4. ) 

X=l.  -( SIGMA*THIK*3.1416*RBUL*RBUL/E0I  ) 

FO=X*EOI 
C  CALCULATION  OF  ENERGY  GIVEN  UP  IN  SHOCK  PHASE 

IF  (1FLAG.GT.0)  GO  TO  3 

X6M=VI/C0 

XE=PIF2(XBM,XXM,32,XX) 

EO=XE-EO 

GO  TO  4 

3  E0=E03 

4  RO=(EO*32.2*144./(RHO0*C0*C0*2.*3.1416) )  ** 0.3333 3 
C 

C      INITIAL  VALUE  OF  ENERGY/VOL. 
C 

ETO=2.*EO/(RO**3.0) 
C 

C  INTEGRATE    ENERGY     EQUATION    TO    FIND    SHOCK    MACH    #, 

C  DENSITY    AND    PRESSURE    RATIO. 

C 

CALL    ENERG    ( IFL AG , ETO ,C0, RHOO, XN , VMS , PS  I  10 , F 10) 

VSO=VMS 

VSI=0.01 

CALL  DEVF1  ( VSO , VSI , X N, DODMO ) 
C 

C      COMPUTE  INITIAL  TIME  FOR  DISTANCE  RO 
C 

T0=R0/(C0*VS0*12.0) 

c 

C      ASSIGN  SUBSEQUENT  VALUES  OF  SHOCK  MACH  NO.  AT  WHICH 

C      SHOCK  RADIUS  WILL  BE  COMPUTED 

C 

DO  10   1=1, N 

IF  (VSO-VS(D)   1,2,2 

1  GO  TO  10 

2  M=N-I+1 
XMS( 1)=VS0 
GO  TO  11 

10  CONTINUE 
C 

C      COMPUTE  SUBSEQUENT  VALUES  OF  SHOCK  RADIUS 
C 

11  RS(1)=R0 

DO  20   J=2,M 

XMS( J)=VS ( I) 

VSE=XMS( J) 

CALL  SUBEN  ( I FL AG , CO , RHOO , XN, VSE , ET ) 

RS( J)=( 2.*E0/(ET) )**0.3333 

1  =  1  +  1 
20  CONTINUE 

K  =  M-1 
C 

C      COMPUTE  DELTA  TIME  CORRESPONDING  TO  SHOCK  MACH  NO. 
C      AND  SHOCK  RADIOS 


C 


DO  30   1  =  1  ,K 
DRS( I)=  RSII+1)-RS( I) 
ARMS(  I )=2./(XMS(  I  )+XMS(  If  1 )  ) 
30  DT( I )=DRS( I )*ARMS( I) /(CO* 12.0) 
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T(1)=T0 

DO       40       J=2,M 
40    T( J)=    T( J-1)+DT( J-l  ) 

IF(RHO0.LT.50.3)     WRITE(6,104) 

IFCRH00.GT.61.0)     WRITE(6f102) 

WRITE  (6,103)   EO 

WRITE  (6,110) 
C 

C      PRINT  SHOCK  MACH  NO.  AND  SHOCK  RADIUS  V S  TIME 
C 

DO  50   1=1 ,M 
50  WRITE  (6,111)   T ( I ) ,XMS( I ) ,RS( I ) 

100  FORMAT   (8E10.0) 

101  FORMAT   (3E10.0) 

10?  FORMAT  (1H,// ,?9X, 'ENERGY' ,15X, • FLUID  WATER') 

103  FORMAT  (29X, 1E20.7) 

104  FORMAT* 1H,//,39X, « ENERGY' ,15X, 'FLUID  FJEl.') 

110  FORMAT  (1H,//,39X  ' T I  ME ' , 1 7X , ' MACH •  ,  1 7 X  ,  • RADI US '  ) 

111  FORMAT   (29X,3E20.7) 

200  FORMAT  ( //  10E 1 2 . 5/ / ) 

201  FORMAT  (15X,2E20.7) 

202  FORMAT*  15X,1E20.7) 

203  FORMAT  (I2,2E10.0) 

204  FORMAT  ( // 13X I  2 , 20XE 1 2 . 5// ) 
C 

C 

C      PART  II 

C 

C 

C  INITIALIZE    VARIABLES    AND    ASSIGN    VALUES    USED    WITH 

C  NON-DIMENSIONAL    PRESSURE    DISTRIBUTION    FUNCTION 

C 


BBB  =  T(M)*1.E06H 

.05 

L  =  BBB 

TT(1)=0.0 

XM(1)=XMS(1) 

R(1)=R0 

F2(l )=F10 

PSI2(1)=PSI10 

D0DM2( 1 )=DuOMO 

Q(1)=3.*(PSI10-1 

.0) 

B=-0.25 

I=1.0/DC  +1.05 

C( 1)=0.0 

WRITE  (6,301) 

WRITE  ( 6,302 )  TT(  1  ) 

WRITE  (6,303) 

DO  60    J=2,M 

C (J)=C( J-l )+DC 

60 
C 

C      COMPUTE  N0N-D.1MENSIQNAL  AND  DIMENSIONAL  PRESSURE  FOR 
C      INITIAL  SHOCK  RADIUS. 
C 

DO  70  J  =  l,  I 

F(lt J)-F2(  1)+PSI2(1)*(C( J)**(Q(l)+2.0)-1.0)*(F2( 1)/ 
1PSI2(1)-B*(F2(  1  )  +XM(  1  )*D0DM2(  L)))/(Q(  1)1-2.0) 

P(1,J)=F(1,J)*RHG0*C0*C0*XM( l)*XM(l)/(  32 .2*144.  ) 
70  CONTINUE 

K=l 

JJ=1 

SUM=0.0 

CALL    FORCE(FOR,C, P,R, I ,JJ,K) 

DO    501    J=l, I 

SUM=SUM+FOR(  J) 

WRITE  (6,304)     C  (  J  )  ,  P(  1  ,  J  )  ,  F  (  1 ,  J  )  ,  FOR  (  J  ) 
501    CONTINUE 

TOFOR( 1 )=SUM 
C 

C      ITERATION  TO  SOLVE  PRESSURE,  SHOCK  RADIUS,  SHOCK  MACH 
C      NO.  VS  TIME 
C 

DO  80   K=2,L 
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TT(K)=TT(K-1)+DELT 


21,21,22 


C 

c 
c 

c 


IF  (TT(K)-TO) 
21  XM(K)=XM(1) 
R(K)=R( 1 ) 
F2(K)=F2( 1  ) 
PSI2(K) =PSI2(1 ) 
DODM2(K)=DODM2( 1 ) 
Q(K)  =  Q(l ) 
DO  55   J=l 
P(2,  J)  =  P(  1 
GO  TO  80 
X=TT(K) 


INTERPOLATION  OF  INPJT  DATA  AND 
VARIABLES  IN  PRESSURE  FUNCTION 


55 
22 


.1 

t  J) 


PART  I  DATA  FOR 


90 


XM(K) =PIF2(X,T,M,XMS) 
R(K)=PIF2(X,T,M,RS) 
VSS=XM( K) 

CALL  SHOCK  ( XN, VSS , PS  I  3 , F 12 ) 
PSI2(KJ =PSI3 
F2(K)=F12 

CALL  DEVF1  ( VSS , VS I ,XN ,D0DM3 J 
D0DM2(K )=D0DM3 
Q(K)=3.*(PSI2(K)-1.0) 
WRITE   (6,301) 

WRITE   (6,302)  TT ( K) , XM( K) , R ( K ) 
WRITE   (6,303) 
DO  90   J=l ,1 

F(2,J)  =  F2(K)+PSI2(K)MC(J)**(Q(K)+2.0)-l.G)*(F2(K)/ 
1PSI2(K)-B*(F2(K) +XM(< )*D0DM2(K) ))/(Q(K)+2.0) 

2*144.  ) 


:C0*XM(K)*XM(K)/(32 


P,R,  I  ,  JJ,K) 


C 
C 
C 


502 

80 


P(  2,J)=F(  2,J)*^H00*C0: 

CONTINUE 

JJ  =  2 

SUM =0.0 

CALL  FORCE(FOR,C 

DO  502  J=l ,  I 

SUM=SUM+FOR( J) 

WRITE (6, 304}  C(J),P(2,J),F(2,J),FOR(J) 

CONTINUE 

TOFOR(K)=SUM 

CONTINUE 

ACOUSTIC  EXTENSION  OF  THE  MACH  WAVE. 


T  T  (  K  )  ,  X  V  (  K )  ,  R  (  K  ) 


*RH00*CD*C0/(32.2* 
C( J ) ,P(3, J) ,F(2,J 


KK=TMAX+1.05 

LL=L+1 

DO    610    K=LL,KK 

TT(K)=TT(K-1 JOELT 

XM(K)=1 .00 

R(K)=R(K-1 )+C0*DELT*12 

WRITE     (6,301 ) 

WRITE    (6,302) 

WRITE     (6,303) 

DO    600    J=1,I 

P(3, J)=F(2, J) 
600  WRITE  (6,306) 
610    CONTINUE 

WRITE(6,305) 

CALL    CONVER(R,XM,TGFu*tTTtL) 

91  FORMAT     (3E15.5) 

92  FORMAT     (1E15.5) 

301  FORMAT     ( 1H,//  ,36X 

302  FORMAT     (29X,3E2D. 

303  FORMAT(  1H,//,39X, 

304  FORMAT(29X,4E20.7 

305  FORMAT( 1H,//,36X, 
10F0RCE'  ) 

306  FORMAT     (29X,3E20 
END 


144.*R(<) > 
) 


1  TI  ME  •  ,17X, 'MACH'  ,17X, 'RADIUS'     ) 


7) 

'^/ 

) 

•TIME 

7) 


S',17X,'P',17X,'F',19X,  'FORCE') 
•,17X,'VELO',17X,«  RADUISe , 1 7X, • T 
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FUNCTION  PIF2  ( X , X L I  ST , N, FL I  ST ) 
C 

C      SECOND  ORDER  INTERPOLATION 
C 

DIMENSION  XLIST  (100),  FLIST  (100) 

BLIF  (P,Q,R,S,T)  =  ( ( G-P)*(S-T)/(R-Q)fS) 

IF  (X-XLIST(N) )  2,1,1 

1  I  =  N-l 
GO  TO  5 

2  IFU-XLISTd  )  )  4,4,6 

4  1  =  1 

5  K  =  1 

GO  TO  30 

6  K  =  2 

7  DO  8  I  =  1 ,N 

IF  (X-XLIST( I ) )  9,9,3 

8  CONTINUE 
I=N 

9  1=1-1 

30  BLIF1  =  BLIF(X,XLIST(I ) ,XLIST( 1+1) ,FLIST( I ) ,FLIST( 1+1) 
1) 

10  IF  (K-l  )  11, 11,12 

11  PIF2  =  BLIF1 
RETURN 

12  IF((I+2)-N)  13,13,16 

13  IF  ( (  1-1  )-l)  15,  14,14 

14  IF(A8S(XLIST(I-l)-X)-ASS(XLlST(I+2)-X))16,15,15 

15  L  =  1+2 
GO  TO  17 

16  L  =  1-1 

17  BLIF2  =  BLIF  ( X , XL  I  S T ( I ) , XL  I  ST  (  L )  ,FLI ST (  I  )  , FLI ST ( L ) ) 
PIF2  =  BLIF  ( X,XLIST( 1+1 ) ,XLIST(L J ,BLIF1 ,BLIF2) 

18  RETURN 
END 

SUBROUTINE  FORCE  (  FOR,  C  ,  P  ,<? , I,JJ,K) 
C 

C  THIS  SUBROUTINE  COMPUTES  THE  FORCE  ON  THE  ENTRY  WALL 

C  DUE  TO  SHOCK  PHASE  FOP  A  FRONT  WALL  WITH  A 

C  PRE-PUNCHED  HOLE. 

C  RHRS=  NON-DIMENSIONAL  ^A^IUS  OF  PRE-PUNCHED  HOLE 

C  S=  SLOPE  OF  STRAIGHT  LIME  BETWEEN  PRESSURE  DATA 
C  POINTS 

C  RX=  THE  INTERCEPT  OF  R/RS, NON-D IMENS IONAL  RADIUS, 
C  WHERE  PRESSURE  IS  ZERD 

C  FOR(N)=  F0RCE(L3)  F3^  A\  AMNULUS  WHOSE  OUTSIDE  RADIUS 
C  IS  C(N)  AND  INSIDE  RADIUS  IS  C(N-l) 

C 

c 

DIMENSION  FOR (20) ,C(100),P(2,20) ,R(300) 
RHRS=.375/R(K) 
IF(RHRS-1.)  100,25,25 
100  IF( JJ-1 )  1,1,2 
1  DO  10  N=2, I 

S=(P( 1,N)-P( 1,N-1) )*4. 
IF(P(1,N-D)  3,4,4 

3  IF(P( 1,N) )  5,5,6 

5  FOR(N)=0.0 
GO  TO  10 

6  RX=(S*C(N-1)-P(  1 ,N-1)  )/S 
IF(RHRS.GT.C(M) )  GO  TO  43 
IF(RX.GT.RHRS)  GO  TO  63 

PH  =  S*( RHRS-C(N-1 )  ) +P(  1  ,N-l ) 

FGR(N)=(.5*(P(1  ,N)-PH)*PH)*3.141-6*(C(N)*C(N)-RHRS 
1*RHRS) *R(K)*R(K) 
GO  TO  10 
60  FCR(N)=P( 1 ,N)*.5*3.1416*(C(N)*C(N)-RX*RX)*R(K)*R(K) 
GO  TO  10 

4  IF(RHRS-C(N-1) )  40,40,41 

40  F0R(N)=(.5*(P(1,N)-P(  1 , N- 1  )  ) +P ( 1  ,  N- 1 )  ) *  3 . 1 4 16* ( C ( N ) 
1 *C ( N)-C (N- I ) *C ( N-l ) ) *R( K) *R( K) 
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43 

10 


13 
15 

16 


GC    TO    10 

41  IF(RHRS-C(N) )     ^2,43,43 

42  PH=S*( RHRS-C(N-1 I  )+P<  1 , N- 1 ) 

FOR(M)  =  (  .5*(  P(  U  N)-PH)*PH)*3.1416*<C(N)*C(N)-RHRS 
l*RHRS)*R(K)*R(K) 

GO    TO    10 

FOR (N) =0.0 

CONTINUE 

RETURN 

DO    20    N=2, I 

S=(P( 2»N)-P( 2.M-1) )*4. 

IF(P<  2,N-1)  )     13, 14,14 

IF( P(2,N)  )     15,15  ,  16 

FOR (N) =0.0 

GO    TO    20 

RX=(S*C(N-1)-P(2,N-1 ) )/S 

IF(RHRS.GT.C(M) )     GO    TO    143 

IF(RX.GT.RHRS)     GO    T3    loO 

PH=S*( KHRS-C(N-I ) )+P{ 2,N-l ) 

F0R(N)=(.5*( P(2,N)-PH)+PH)*3.1416*(C<N)*C(N)-RHRS 
i*RHRS)*R(K)*R(K) 

GO    TO    20 
160    FOR(N)=P( 2,N)*.5*3.1416*(C(N)*C(N)-RX*RX )*R(K)*R(K) 

GO    TO    20 
14    I F(RHRS-C(N-1) )     140,140,141 

140  F0R(N)=(.5*(P{2,N)-3(?,N-1))+P<2,N-1))*3.1416*(C(N) 
1*C(N)-C(N-1)*C(N-1) )*R(K)*R(K) 

GO    TO    20 

141  IF(RHRS-C(N) )     142,1^3.143 

142  PH=S*( KHRS-C1N-1 ) ) +P(2,N-1 ) 

FOR  <  N )=(. 5 * ( P ( 2  f  N ) - P H  )  +PH ) *3 . 1 41 6* ( C ( N ) *C ( N ) -RHR S 
i*RHRS)*R(K)*R(K) 
GO  TO  20  \ 

143  FOR(N)=0.0 
20  CONTINUE 
25  RETURN 

ENO 


C 

c 
c 
c 


SUBROUTINE  COM VER ( R , XM , TO FOR , TT , L ) 

THIS  SUBROUTINE  CONVERTS  VELOCITY  AMD  RADIUS 
TO  KM/ SEC  AND  CM  RESPECTIVELY 

DIMENSION  R(300) ,XM(300) , TO FOR < 300) , COMR S ( 300 ) , SHOC V ( 
1300),TT(300) 

DO  510  11=1, L 

CONRS(  II)=R( II)*2.54 

SHOCV( II)=XM(II)*4.9/3.281 

WRITE (6,520)  TT(II) ,SHOCV( 1  I ) ,CONRS( 1 1 )  , TOFORt 1 1 ) 
520  FGPMAT(29X,^E20.7) 
510  CONTINUE 

RETURN 

END 


C 
C 


SUBROUTINE    SHOCK ( XN, VS, PS  I  1 , F I ) 

THIS    SUBROUTINE    C0MPJTE6    THE    NON-DIMENSIONAL     PRESSURE 
AMD    DENSITY    AT    THE    SHOCK    FRONT. 

F(X)=-X+1/(1-B*(X**XN-1)) 

X1=(XN*( VS**2. ) ) **( i/XNi-0.05 
6=1/ (XN*(VS**2. ) ) 

5    DO    10    I =1,200 

Xl=Xl-0.005 

IF(X1  .LT.O.O)     30    TO    20 

Y=F(X1 ) 

IF(Y.LE. 0.001)     CO    TO    40 
10    CONTINUE 
20    DO    30    J=l, 200 

Xl=Xl+0.0025 
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I F ( XI .GT.O.O)  3C  TO  5 

Y=F(X1  ) 

IF(Y.LE. O.OOl )  30  TO  40 
30  CONTINUE 
40  PSI1=X1 

Fl  =  (  (PS  I 1**XN)-1 . )/<XN*(VS**2. ) ) 

RETURN 

END 

SUBROUTINE    DEVF1 ( VS ,VS I ,XN, DODM ) 
C 

C  THIS    SUBROUTINE    COMPUTES    THE    DERIVATIVE     OF    THE 

C  NON-DIMENSIONAL     PARTICLE    VELOCITY    WITH    RESPECT    TO    THE 

C      SHOCK  MACH  NUMBER. 

C 

VSA=VS+VSI 

VSB=VS-VSI 

CALL    SHCCK(XN,\/SA,  PSIA,FA) 

CALL    SHCCK(XN,VSB,PSI8,FB) 

D0DM=(FA-FB)/(2.*VSI ) 

RETURN 

END 

SUBROUTINE    ENER3    ( IFLAG,ETO ,C0 , RHOO,  XN,  VMS , PS  1 10, F10) 
C 

C  THIS    SUBROUTINE    COMPUTES    THE    INITIAL    SHOCK    MACH    NUMBER 

C  BY    INTEGRATION    3 F    THE    EMERGY    EQUATION.       THE    EXTERNAL 

C  SUBROUTINE    QG10     IS    A    GUASSIAN    10    POINT    QUADRATURE. 


C 


EXTERNAL  FCT 

COMMON  BQ 

C  =  2.0 

IF(IFLAG.EQ.O)  C=1.0 

VSG=4.0 

CALL  SHOCK  (  XM  ,  V  S3 , ^ S I  10, F 1 0 ) 

A=3.1416*(RHGO/3  2.2)*(CO*CO/144.0)*(VS3**2.0)*PSI  10 
1*F10*F10*2.0 

BQ=3.*PSI10-1. 

CALL  QG10  (0.0,1 .0, FCT, EP) 

ETG=A*EP*C 

IF(ETO-ETG)  40,100,50 
40  VSG=VSG-.l 

CALL    SHOCK     ( XN , V SG , PS  I  1 0 , F 1 0 ) 

BQ=3.*PSI10-1. 

CALL    QG10     (0.0,1 .0, FCT, EP) 

ETG=A*EP*C 

IF(ETG-ETO)    60,100,40 
60    VSG=VSG+.001 

CALL    SHOCK    ( XN, VSG ,PS 1 10, F 10) 

A=3.1416*(RHO0/3  2.2)*(C0*C0/144.O)*( VSG**2.0)*PSI 10 
1*F10*F10*2.0 

BQ=3.*PSI10-1. 

CALL    QG10     (0. 0,1.0, FCT, EP) 

ETG=A*EP*C 

IF(ETG-ETO)    60,100,130 
50    VSG=VSG+.l 

CALL    SHOCK    ( XN , V  S3 , P S I  10 , F 1 0 ) 

A=3.1416*(RHO0/3  2.2)*(CO*C0/144.0)*( VS3*#2. 0) *PS I  10 
1*F10*F 10*2.0 

BQ=3.*PSI10-1. 

CALL  QG10  (0.0,1 .0, FCT, EP) 

ETG=A*EP*C 

IF(ETG-ETO)  50,100,90 
90  VSG=VSG-.001 

CALL  SHOCK  ( XN , V SG , PS  I  1 0, F 1 0 ) 

A=3. 141 6* (RHO0/ 3 2 .2)-~<CO*CO/144.0)*(VS3**2.0)*PSI10 
1*F10*F 10*2.0 

BQ=3.*PSI10-1. 

CALL  QG10  (0.0,1 .0, FCT, EP) 

ETG=A*EP*C 
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IF(ETG-ETO)  100,100,90 
100  VMS=VSG 
RETURN 
END 

SUBROUTINE  SUBEN ( I FL AG , CO , RHOO , XN , VSE , ET ) 
C 

C      THIS  SUBROUTINE  COMPUTES  THE  ENERGY  BEHIND  THE  SHOCK 
C      FRONT  FOR  A  GIVEN  V.ACH  NO.   THE  EXTERNAL  INTEGRATION 
C      SUBROUTINE  0G1D  IS  A  GAUSSIAN  10  POINT  QUADRATURE. 
C 

EXTERNAL  FCT 

COMMON  BQ 

C  =  2.0 

IF( IFLAG.EQ.O)  C=1.0 

CALL  SHOCK  ( XN, VSE  ,  PS  I E , FE ) 

A= 3. 141 6* (RHOO/ 3 2.  2  )*(CO*CO  / 144.0  )*( VSE**  2.0)  *PS  IE 
1*FE*FE*2.0 

BQ=3.*PSIE-1 . 

CALL  QG10  (0.0, 1.0, FCT, EE) 

ET=A*EE*C 

RETURN 

END 

FUNCTION  FCT(X) 
C 

C      THIS  FUNCTION  SJ3PR03RAM  GENERATES  FUNCTIONS  FOR  THE 
C      INTEGRATION  ROUTINES. 


C 


COMMON  BQ 

FCT=( l.+X**2. )*X**3Q 

RETURN 

END 
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